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Abstract 

Provided that Supcrsymmctry (SUSY) is realized, the Large Hadron CoUider (LHC) 
and the future International Linear Collider (ILC) may provide a wealth of precise data 
from SUSY processes. An important task will be to extract the Lagrangian parameters. 
On this basis the goal is to uncover the underlying symmetry breaking mechanism from the 
measured observables. In order to determine the SUSY parameters, the program Fittino 
has been developed. It uses an iterative fitting technique and a Simulated Annealing 
algorithm to determine the SUSY parameters directly from the observables without any 
a priori knowledge of the parameters, using all available loop-corrections to masses and 
couplings. Simulated Annealing is implemented as a stable and efficient method for finding 
the optimal parameter values. The theoretical predictions can be provided from any 
program with SUSY Les Houches Accord interface. As fit result, a set of parameters 
including the full error matrix and two-dimensional uncertainty contours are obtained. 
Pull distributions can automatically be created and allow an independent cross-check of the 
fit results and possible systematic shifts in the parameter determination. A determination 
of the importance of the individual observables for the measurement of each parameter 
can be performed after the fit. A fiexible user interface is implemented, allowing a wide 
range of different types of observables and a wide range of parameters to be used. 



1 Introduction 



Fittino [1] is a program implemented in C++, which extracts the parameters of the SUSY 
Lagrangian from (simulated) measurements at future colliders such as the LHC and the ILC 
in a global fit. It has been created in the context of the Supersymmetry Parameter Analysis 
(SPA) [2, 3] project. In Fittino, no a priori knowledge of the parameters is assumed. In 
contrast to the approach pursued in [4, 5], tree-level relations among observables and subsets 
of SUSY parameters are used to obtain start values for the x^"fit- Fits in parts of the parameter 
space, using only observables directly depending on the fitted subset of parameters, are used to 
refine the tree-level parameter estimates. Alternatively, a Simulated Annealing algorithm can 
be used to globally minimise the for all parameters with respect to all observables. Finally, 
a global fit and an uncertainty analysis provide the best parameter values, their uncertainties 
and their correlations. Additionally, two-dimensional graphical contours of the parameter 
uncertainties can be obtained, pull distributions and distributions can be automatically 
created and an analysis of the importance of individual observables for the determination of 
each parameter can be performed. 

Fittino allows to perform a simultaneous fit of any combination of 24 parameters of 
the low energy MSSM Lagrangian and 7 Standard Model (SM) parameters. Alternatively 
high-scale mSUGRA, GMSB or AMSB parameters can be fitted. As observables, measure- 
ments of present or future collider experiments can be used in the fit, such as masses, cross- 
sections, branching fractions, widths, products of cross-sections and branching fractions, ra- 
tios of branching fractions and edges in mass spectra. Correlations among observables can be 
specified. This allows a realistic test of the precision of the parameter determination of the 
low-energy MSSM theory at future colliders. The data of measured and simulated observables 
and their uncertainties (if desired including theoretical uncertainties) and correlations is en- 
tered by the user. The theory prediction for these observables is calculated by a theory code, 
which is interfaced to Fittino via the SUSY Les Houches Accord (SLHA) [6]. In the present 
implementation, SPheno [7] is used as an example to calculate the observable predictions from 
the set of parameters. 

The user interface of Fittino allows the user to specify any number of observables and 
parameters as given in the list above. The behaviour of Fittino can be steered with a wide 
range of options. Start values for the global fit can either be determined by Fittino, using 
tree-level estimates, several sequential fits in parts of the parameter space (subsector fit) or 
Simulated Annealing, or they can be specified by the user if desired. 

Fittino has been tested for different MSSM scenarios. The fit strategies of Fittino, both 
with subsector fits and with the Simulated Annealing algorithm, provide stable parameter 
reconstruction using simulated measurements from ILC and LHC. The results of these fits are 
reported in [8]. 

This paper is organised as follows: First the general concepts of MSSM parameter determi- 
nation are introduced, followed by a description of the fit program Fittino. A documentation 
of Fittino is given starting from Section 3, where the syntax and the commands of the input file 
are discussed. In Section 4, a short guideline for how to start a fit with Fittino is given. This 
is followed by a description of the output in the remaining chapters. In Section 6, a summary 
and a short description of the results obtained with Fittino for various MSSM scenarios [8] 
are given. All information provided in this article refers to Fittino version 1.1.1. 
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2 The Fit Program Fittino 



In the presence of experimental uncertainties which are much larger than loop effects in a 
theory, parameters can be measured using analytical tree-level relations among parameters 
and observables. However, for very good experimental accuracy this method is depreciated for 
a correct parameter and uncertainty determination, since parameters tend to be systematically 
off their true values and since correlations can not be fully taken into account. This plays a 
role if the relative uncertainties of the measured observables Oi are smaller than the largest 
contribution from loop effects 

^Qmeas AO'°°^ 

as in case of the future measurements of the ILC. Therefore in this case all available loop 
corrections have to be taken into account, in order to achieve the highest possible precision. 
This means that all observables Oi are treated as functions of all parameters Pj: 

Observable Oi = / (all parameters Pj). 

Additionally, in order to account for the uncertainties from the limited precision of SM param- 
eters (parametric uncertainties), the SM parameters can be fitted simultaneously with MSSM 
parameters. This approach also allows to extract the full correlation among all parameters. 
No bias is introduced due to a priori assumptions on fixed parameters. 



2.1 The Fittino Approach 

The aim of Fittino is the unbiased determination of the parameters of the MSSM Lagrangian 
-Cmssm, obeying the following principles: 

• No a priori knowledge of SUSY parameters is assumed (but can be used if desired by 
the user). 

• All measurements from future colliders could be used. 

• All correlations among parameters and all influences of loop-induced effects are taken 
into account. 

It is impossible to determine all 105 possible parameters of >Cmssm simultaneously. There- 
fore, assumptions on the structure of -Cmssm are made. All complex phases are set to 0, no 
mixing between generations is assumed and the mixing within the flrst two generations is set 
to 0. Thus the number of free parameters is reduced to 24 (MSSM-24). Further assumptions 
can be specified by the user. Observables used in the fit can be 

• Masses, limits on masses of unobserved particles 

• Widths 

• Cross-sections (momentarily in e"'"e~ collisions only) 

• Branching fractions 

• Edges in mass spectra of decay products. For example in slepton decays — > 
//"'"XiM~/ti) the lepton energy spectrum has a box-like shape. Its edges are correlated to 
the masses of the ft and the Xi- By using the edge positions instead of the reconstructed 
masses, correlations among the observables can be reduced or omitted (see e.g. [9]). 
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• Products of cross-sections and branching fractions. For most processes, it is hard to 
measure total cross-sections. Therefore an observable can be formed as a product of a 
cross-section and several branching fractions. 

• Ratios of branching fractions. If the total width of a decay can not be measured, the 
individual measured widths can be all normalised to one width and thus form ratios of 
branching fractions. 

Correlations among observables and both experimental and theoretical errors can be speci- 
fied. Theoretical uncertainties can be important, if they are in the order of magnitude of the 
experimental uncertainties. Both SM and MSSM observables can be used in the fit. Para- 
metric uncertainties of SUSY observables can be taken into account by fitting the relevant 
SM parameters simultaneously with the MSSM-24, the mSUGRA, the GMSB or the AMSB 
parameters. For the interface between Fittino and the code providing the theoretical pre- 
dictions, the SUSY Les Houches Accord [6] (SLHA), is used. The SLHA is a format for a 
text-file based interface between spectrum calculators, event generators and other programs in 
the context of the MSSM. Any theoretical code compliant with SLHA can be easily interfaced 
with Fittino. In the current implementation, the prediction of the MSSM observables for a 
given set of parameters is obtained from SPheno [7]. A Simulated Annealing algorithm based 
on [10, 11] and MINUIT [12] is used for the fitting process. 

2.2 The Iterative Fit Procedure 

The full MSSM-24 parameter space in Fittino, consisting of MSSM parameters plus SM pa- 
rameters, cannot be scanned completely, neither in a fit nor in a grid approach. Therefore, in 
order to find the true parameter set P in the presence of nobs measured observables Of^^^ in 
a fit by minimising a function (x^ = Yl^=\^ i^f^^^^ — Of^*^" {P)y / af) , it is essential to begin 
with reasonable start values, allowing for a smooth transition to the true minimum. The three 
different tasks in the parameter determination are: 

1. Tree-level estimate 

Using a small number of observables, the initial values for the parameters are determined 
on tree-level. 

2. Finding the central values of the parameters 

Starting from the initial values of the parameters, the configuration of parameter values 
with the smallest total x^) i-e- the best agreement between theory predictions and 
measurements, has to be found. 

3. Uncertainty determination 

Once the optimal parameter values are found, the evaluation of their correlations and 
uncertainties is performed. 

The iterative procedure of Fittino to determine the start values for the fit, the optimal 
parameter values and the parameter uncertainties is displayed in Fig. 1. In a first step, the 
SUSY parameters are estimated using tree-level relations. Then several subsequent subsector 
fits in parts of the parameter space can be used to improve the tree level estimates before the 
global fits. This is called the subsector fit method. Alternatively, a Simulated Annealing [10, 
11] algorithm can be used to find the minimum. In a third step a global fit of all parameters to 
all observables is used to fine-tune the optimal parameter values and determine the parameter 
uncertainties and correlations. 
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fittino.in: Observables O, Uncertainties AO 



V 



Tree-level estimate: P = Tree_level_function (O) 




Tree-level 
estimate 



Simulated Annealing: 
X minimization with all 
parameters P and all observables O 



Simulated 
Annealing 



I 




Fit parameter m^ 
to observable mA 



Fit tanP and (j, 
to gaugino 
observables 



Fit tanp, slepton 
mass par. and to 
slepton observables 



Scan Xjj and Xt 
using squark 
observables 



Fit tanP, squark 
mass par. and Xj \, to 
squark observables 

Fit tanP and Xj ^ 
to all observables 



UseSimAnnWhile : Shake fitted 
parameters with Simulated Annealing 
I 

Global fit of all parameters P to 
all observables O, using the start 
values determined by the previous 
steps 



Loop: NiimberMinimizations 



Re-evaluate the uncertainty matrix 
using the HESSE command 



Refine the uncertainty matrix 
and determine asymm. eiTors 
using the MINOS command 



Determine 2-dimensional 
uncertainty contours 



Store the fit result in fittino.out 



Separate 
in parts of^ 
the W 

parameter 
Space using 
restricted 
sets of 
observable 



I 




Global fit 




Uncertainty 
calculation 




Figure 1: Flow chart of the iterative fit procedure of Fittino. Output and input files are shown 
in light courier, boolean operations in bold courier and explanations are shown in times. 
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2.2.1 Tree level estimates for the parameters 

In order to get the start values for both the subsector fit method and the Simulated Annealing, 
tree-level relations of the form 

Observable = /(Parameters) 

are inverted to the form 

Parameter = /'(Observables) 

and used to estimate the parameters. This is done in the gaugino, the slepton and the squark 
sectors. 

1. /u, myi, tan/3, Ml, M2, M3 are determined from the gaugino and Higgs sector using for- 
mulae from [13, 14]. In order to extract these parameters, information from chargino 
cross-sections is needed, which enters in form of the chargino mixing angles cos2i;^j;, and 
cos 2(j)ji. These pseudo-observables are just used for the determination of the start val- 
ues, no use is made of them for the fit. They can be approximately determined from 
chargino cross-sections at the ILC at different polarisations. If kinematically accessible, 
mApoie can be directly measured. The start value of the parameter m-Arun is set to 
mApoic ill this case. For the other parameters, the following individual calculations are 
performed. For the full form of the relations, see e. g. [7]. Based on the diagonalisation 
of the chargino mixing matrix the mixing angles <j)L and can then be used to derive 
the tree-level estimates of tan/? and M2: 

= mvK (5^ + A(cos 2<^L + cos 2<^/j)) 2 (1) 

_ /^ l + A(cos2(/)fl-cos20z.) ^^ 
^""""^ ~ Vl-A(cos2</>^-cos20L)y' 

M2 = mvF (i; - A(cos2(/)L + cos2(/)r))5 (3) 



sign(//) = g 2 A^l I ■ o/^ ^) 

M3 = mg (5) 



usmg 



X2 Xi 

2m^ 

2 2 
m ± — m ± 



(6) 



A = , (7) 

After these parameters have been determined at tree-level, a tree-level estimate of Mi can 
be calculated from the neutralino mass matrix eigenvalues. They depend on ^, tan/3, M2 
and Ml. The first three of these parameters have been determined already in the first 
step. Thus the neutralino system can be used to finally determine Mi. In fact, the neu- 
tralino system is the only sector where Mi can be measured directly. The characteristic 
equation of the neutralino mass matrix then is used to determine Mi [14]. 

2. Xt, Xb, Mq, M[/, M/) are determined from the squark sector masses, using formulae 
from [15]. The trilinear coupling parameters A are set to zero for the determination 
of the squark mass parameters. Xt = At — fi/tan (3 and X]~, = A\j — ^tan/3 are chosen as 
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fit parameters instead of A because their correlation with tan/5 is reduced. After tan/3 
and fi have been determined in the previous step, the following tree-level relations can 
be used. As an example, only the third generation is shown explicitly. 



2 2.2«s o 1 



= -mz' cos2p{- -- sin' ew)-mt + -{ml + ml) (8) 

= mz^cos2(3{^-^sm^9w)-ml + ^iml^+mlj (9) 

Mf ^ = -mz^ cos 2/?^ sin^ Ow + \ ("^t'l + ) (^0) 

M^^ = mz2cos2/3^sin2 0H/-rnb + ^(Ki +"^62) (1^) 

Xt = -/i/tan/3 (12) 

Xb = — ^tan/3 (13) 



The parameters of the other generations can be obtained analogously. The initial deter- 
mination of and X\, is only very rough and therefore refined in an additional step, 
which is explained below. 

3. Xr, Ml, Me are determined from the slepton sector masses, using formulae from [15]. 
Their calculation is analogous to the calculation of the squark parameters in the previous 
item. The determination of X^ is refined later as explained below. 

In order to estimate the uncertainty of this determination of the parameters from tree-level 
formulae, the calculation is repeated 10 000 times with observables randomly smeared within 
their uncertainties according to a Gaussian probability distribution. The starting value of the 
fit is the mean of the distribution of each parameter. From the variance of the calculated 
parameter distribution the initial uncertainty is estimated. 



2.2.2 Finding the optimal parameter values 

A global fit of a large number of MSSM parameters with starting values derived from 
tree-level relations will most likely fail to find the global minimum, since the parameter 
space contains many local minima. Therefore the tree-level estimates have to be improved 
before the global fit. This can be done by using knowledge about the dependencies of sets of 
individual parameters on sets of individual observables. This is done in the fast subsector fit 
method. However, also in this case some local minima can be too deep to allow MINUIT to 
escape from them. Therefore the slower but very stable Simulated Annealing algorithm has 
been implemented alternatively. 



Subsector Fit Method Since the possibility of mixing in the third generation has been 
neglected in the calculation of the tree-level estimates, the parameters Xt, Xy,, X^ are only 
roughly initialised. A global fit with these starting values would most likely not converge. 
Therefore next the estimates from the slepton sector are improved by fitting only the slepton 
parameters X^-, Ml, Me to the observables from the slepton sector, i. e. slepton masses, 
widths and cross-sections. Observables not directly related to the slepton sector can degrade 
the fit result, since parameters of other sectors are likely to be still wrong. In such a case a 
parameter of the slepton sector will be pulled into a wrong direction, in order to compensate 
for the wrong parameters of other sectors. All parameters not from the slepton sector are fixed 
to their estimated tree-level values. In this fit with reduced number of dimensions MINUIT 
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Vary each parameter at a time 
If x^improves: accept new parameter 
If degrades: accept new point if 
allowed by Boltzmann-statistics with 
temperature t 

else discard parameter change 



J 



I 



Vary each parameter 20 times 
I 



Re-adjust the parameter step sizes 



J 



Calculate new temperature: t = r t 
Reset parameters to point with best y} 

I 



Figure 2: Structure of the Simulated Annealing algorithm. The four different loops are indi- 
cated by the grey numbers. 

can handle the correlations among the parameter better than in a global fit with all parameters 
free. 

Then the third generation squark parameters are improved by only fitting 
Xt, X];,, Mq, Mu, Md to the observables of the squark sector, masses, widths and cross- 
sections. All other parameters are fixed to their previous values. 

After this step the correlations among tan (3 and the third generation slepton and squark 
parameters are still not optimally modelled. Therefore another intermediate step is introduced, 
where tan (3, Xt, X-^, Xr and ^ are fitted to all observables and all other parameters are 
fixed to their present values. 

The subsector fit method does not yield a parameter set which corresponds to the global 
minimum of the x^- But starting from the result of the subsector fit method, the subsequent 
global fit (see Section 2.2.3) is able to converge to the global minimum. 

Simulated Annealing In the method of Simulated Annealing, the surface in the MSSM 
parameter space is treated as a potential. A temperature t is defined, which according to 
a Boltzmann-distribution sets the probability of a movement in the parameter space which 
corresponds to an upward movement in the potential, i.e. which results in a larger than the 
parameter setting in the step before. In the process of Simulated Annealing, the temperature 
is reduced step by step. This ensures that large movements over high potential barriers are 
possible in the early stages of the annealing process, in order to escape from local minima. 
Later, the temperature is reduced and the parameter setting is thus forced into areas with low 
X^. The individual steps of the algorithm work as follows: 

First, the initial temperature to is determined by varying each parameter 10 times randomly 
within its estimated uncertainty. The estimated uncertainty is determined from the tree-level 
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estimates. Then half of the standard deviation 



of the values obtained is taken as initial temperature. This choice is large enough to allow 
the algorithm to scan large parts of the parameter space and to escape from local minima. 

The core of the Simulated Annealing algorithm is schematically shown in Fig. 2. In the 
innermost loop (loop No. 1 in Fig. 2), each parameter Pj is varied at a time by a random 
variation within the estimated uncertainties of each parameter APj. These uncertainties are 
initially taken from the tree-level estimates and later refined for each parameter individually, 
as described below. After the parameter variation, the the new parameter setting is 

evaluated. If the new value Xi lies below xf-n where Xi-i is the x'^ value of the last accepted 
parameter set, the new parameter value is accepted and the next parameter is varied. If the 
new value xf is the best x^ value achieved so far, Xi and the accompanying parameter set are 
stored. If Xi > Xi-n then the new point is accepted on the basis of a Boltzmann distribution: 
If a random number pe [0, 1] fulfils 



2 2 



p < e 

then the new parameter value is accepted despite yielding a worse fit than the previous pa- 
rameter set. This ensures that the algorithm is able to escape from local minima. 

In loop 2 the variation of each parameter at a time is repeated 20 times. On each variation 
the number of accepted (nacc) and rejected (ntot — 'T'acc) parameter choices is counted for each 
parameter. After loop 2, the estimated parameter uncertainty APi is iteratively re-adjusted 
such that nacc/^tot lies in the range of 0.4 to 0.6. This is achieved by the transformation 

/ iw _ 6 

AB, > 0.6 : Ai^'^^* = AP"^'^ f 2^:^i2L_ + i 



AP. < 0.4 : AP 



0.4 

^pold 



^ 4 ]1^£C_ 

O ' "tot _|_ 1 

^ OA ^ 

In loop 3 this re-adjustment of the parameter uncertainties is repeated 3n times for n param- 
eters, but at least 60 times in case n < 20. The re-adjustment of the parameters variation 
ensures that for each temperature t the random parameter variation tests a parameter range 
with values of Xi ~ xf-i ~ ^• 

After loop 3 is completed, the temperature is reduced by the multiplication with a factor r, 
which is constant during the algorithm. At the same time, the parameter values for the next 
iteration are set to the parameter set with the lowest found so far. The whole procedure 
is repeated until one of the two exit criteria are met. Either, the allowed number of calls to 
the theory code is reached, or the variation of the x^ value within loop 3 is smaller than the 
variation within the last four lowest x"^ values found, and the variation of all values is smaller 
than 0.0001. In each case the Simulated Annealing algorithm exits with the parameter set 
which yields the best x^ found so far. 

The Simulated Annealing algorithm has been tested on fits with 20 free parameters in the 
scenarios SPSla' [3] and SPS7 [16]. In each case, the convergence to the true minimum of the 

is good and very stable. However, while being very robust, this algorithm is considerably 
slower than the subsector fit method. In Fig. 3, some distributions illustrating the behaviour 
of the Simulated Annealing in the SPSla' scenario are shown. In Fig. 3 (a), the variation of 
tan P with the number of calls to SPheno is shown. In the beginning, a large range of tan /3 
is scanned. For lower t at a larger number of calls, the parameter variation Atan /? is reduced 
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Figure 3: The parameter determination using Simulated Annealing in a fit with 19 parameters 
in the SPSla' [3] scenario. In (a), the tested parameter values of tan/3 are shown versus the 
number of calls to SPheno. In (b) the same evolution is shown against logi(^{t), where t is the 
temperature of the Simulated Annealing. In (c) and (d), the same projections are shown for 
-'^bottom) which in contrast to tan/3 is only weakly constrained by the measurements. 

along with t and convergence is achieved. Thus Simulated Annealing is very effectively coarsely 
scanning a large parameter space at high t, while it is also finely testing the area close to the 

minimum. The same evolution is shown in (b) against the temperature t. In (c) and (d), 
the same evolution is shown for the bottom-squark mixing parameter Xbottonn which is much 
weaker constrained that tan f3 and therefore takes longer to converge. 

Simulated Annealing is a robust and reliable algorithm, which is better capable of escaping 
from local minima than the subsector fit method. Also it is much less dependent on the starting 
values obtained from the tree-level estimate. However, it is more computing time consuming, 
since a high initial temperature is necessary to overcome the local minima near the tree-level 
estimates of the parameters and find the true minimum. 

2.2.3 Determining the parameter uncertainties and correlations 

After the search for the optimal parameters, all parameters are released and a global fit is 
done, using the method MINIMIZE in MINUIT. In case the subsector fit method has been used 
for the improvement of the tree-level estimates of the parameters, this step is also necessary 
to find the optimal parameters, which are generally not yet perfectly determined after the 
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Table 1: Commands in the Fittino steering Ble. All commands given in this table set the input 
observables and the parameters used by Fittino. The full form of the commands is given in 
the text. The correlations must be specifyed after all observables. The steering flags can be 
found in Tab. 2. 



Command 


Explanation 


Keys before any line 


# 


Comment line 


nof it 


The observable after nof it is not used in the fit 


Observables 


mass<name> 


Mass of the particle <name>, see Tab. 3 


edge 


Position of an edge in a mass spectrum 


sigma 


cross section of a process 


BR 


Branching fraction 


width 


Total width of a particle 


limit 


Limits on particle masses 


LEObs 


Low energy observable 


sin2thetaW 


The value of sin^ 6w 


cos2phiL 


The value of cos l 


cos2phiR 


The value of cos 2(/)/j 


xsbr 


Product of cross-sections and branching fractions 


brratio 


Ratio of two branching fractions 


brsum 


Sum of branching fractions 


Correlations 


correlationCoef f icient 


Correlation coefficient of two observables 


Parameters 


f itParameter 


Name and eventually value of a fitted parameter 


f ixParameter 


Name and value of a fixed parameter 


universality 


Specifies which parameters are unified 



subsector fits. If Simulated Annealing with a sufficient number of iterations is used, then the 
optimal parameters after Simulated Annealing are generally very close to the true parameter 
values. 

If the global fit converges, a MINOS error analysis is performed, yielding symmetrical and 
asymmetrical uncertainties, the full correlation matrix and 2D fit contours. 

3 The Fittino Steering File 

Fittino is controlled using an ASCII input file whose name is passed to Fittino as parameter. 
If none is given, fittino . in is assumed. An example for this file is given in Section A.l. The 
user can specify the observables, their values and uncertainties, their correlations, the fitted 
and fixed parameters and universalities among generations. Flags can be used to control the 
general behaviour of Fittino. A list of the available commands is given in Tables 1 and 2. In 
the following, the syntax of the commands is explained. 
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Table 2: Flags in the Fittino steering file. The commands given here determine the fit strategies 
and outputs of Fittino. The full form of the commands is given in the text. The correlations 
must be specified after all observables. The other commands used for the Fittino input can 
be found in Tab. 1. 



Flags 


Fi tMnHpl 


Model tvne siirh as MSSM rnSTIORA etr 


LoopCorrect ions 


UsG full loop corrections as provided, by SUSY calculcitor 


ISR 


SwitcVi on TSR 


TTctf^nn Trf^n ^1" a T"t" Vp n iif^c; 

WoCVJX V dlkJ U CLX U V CL J. LLC O 


kJ UCll Ij li. '^±11 jJdl dlliC/ uCl V fX± LICO 111 J L Kj± CLJ. CUilC L> C J. 


"Fi-t-AI IVii -rani- 1 \t 


H tI" q 1 T~\ is T'Q Tn P't'PT'G is'l' ("inr'ti 
X IL dli JJcxl clliiC LjCI o dl UliCt; 


V.^CLX^£^ LLJ. J.LJ J.C> Li 


1 1 r*!! 1 p "hp mill ri 1 til'T'i r\i 1 1'lon t: 

V_>cXiv^ LlldUC JJ LLll LllO Li 1 U U. tHJliO 




V^diCLlldUe lliQlVlClUcli i-iX CUllLilUULlOlio 


Rnnnrl QrinY 


Ot- U UU U.ll(J.o Cll J t dliU. -/V |-j 




^r'ssn Xj finn a i nPToyp ttH'itio' 


Qo-nF-i -|-Tan"hy 
ot^ ur X i-> 1 diiUA 


X Ci IL*! Ill oCJJdldLC 111) Ul Udll fJ-f IVlj^^ ^ ivjj-^ 5 -^t dllLl 


OcUr X L i dliUrlU. 


jrciioiiii oepdidije iiu ui Udii /j diiu. tx 


oepr xx-uiii 


JreilOiIil aepdldLe nt OI /'t^Arun ''''Apole 


V^dXC UXd L UX 


optJciiy Cdic uidiui iOi Liieoi y pi cciic iiuiib 


Tic* Ivl T m 


Uoc IVliiN^O CiiOi CdlCUldllOil 


TTc!o"HoC! O O 

u o t^nt; o o 


1 G£i H H ^ H , C^TTCW m Q'l'T'TV" r'Q r'll Q'I'l/^iVl 

oC lill/OOXJ ei 1 Ui lildLi lA. l.^diC UldllUli 


iMumuei uii ixrLxmxzcLT. xoiis 


uiiiuei oi iiiiiiiiiiibciLioii oiepb 
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MaxCallsSimAnn 


Maximum number of calls in Simulated Annealing 


InitTempSimAnn 


Initial temperature of the Simulated Annealing 


Verbose 


Verbose mode 


ScanParameters 


Scan one- or two-dimensional parameter space and write 




X^-surface to ROOT file 


scanParameter 


Specify parameters to scan, the scan range and the 




number of scan steps for ScanParameters 


Perf ormFit 


Controls whether fit is performed or not 


RandomGeneratorSeed 


Seed for the random number generator 


MaxCalculatorTime 


Maximal allowed time for the SUSY calculator code 



Keys The commands in this group act on the contents of the line after the key. Available 
keys are: 

• # 

Comment line. Everything after # is ignored by Fittino. 

• nofit <observable> <value> +- <uncertainty> 
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This key can be used to specify observables which shall only be used in the step of the 
initialisation of parameters using tree-level relations. Typical examples are the chargino 
mixing angles cos2(/)i and cos2(/)/j. 

Observables The commands in this group specify the measured (or simulated) observables 
and their values and uncertainties. The following types of observables are available: 

• mass<name> <value> +- <uncertainty> [ +- <theo_uncertainty> ] 

Specifies the mass of the particle <name>. All available particles are listed in Tab. 3. If 
no uncertainty is given, the particle mass is not used in the of the fit. If more than 
one uncertainty is given, the uncertainties are added in quadrature. This is useful in 
case large theoretical uncertainties on the predicted particle mass exist. For the time 
being it is assumed that treating the theoretical uncertainties as uncorrelated with the 
experimental error and a Gaussian probability distribution is reasonable. 

If no unit is given in <value> and <uncertainty>, they are assumed to be given in GeV. 
Other supported units range from eV to EeV. 

• edge <type> <massl> <mass2> [more masses] <value> +- 
<uncertainties> alias <alias> 

Specifies the edge in a mass spectrum. Since SUSY particles tend to decay in cascade 
decays, the masses of intermediate particles can often be reconstructed from edges in 
mass spectra. If the edges are transformed into masses before the fit is made, the cor- 
relations among the reconstructed masses from one spectrum have to be specified using 
correlationCoef f icient. The command edge offers the simpler and more straight- 
forward possibility to use the edge positions in the mass spectra directly in the fit. 
Momentarily, the following types <type> are available: 

1 nil + 

2 |mi — 

3 (m|-m§)(m§-mf) 

4 (m|-m|)(m§-mf) 

5 [{ml + ml){ml - ml){rnl - mf) 
— {ml — 777,2)1/ {m'2 + m|)^(m| + mf)"^ — 16m'2m^m'f 
+2m|(m| — ■m'^){m\ — m\)]/ {Am'l'm^) [17] 

The list of formulae can be easily extended. Here and in most of the following commands 
<alias> is an integer number which can be used to unambigously identify the input if 
it is used in other commands, such as correlationCoef f icient. 

• sigma ( <initial_state> -> <f inal_state_particles> , <Ecms> , <polarisationl> , 
<polarisation2>) [ <value> +- <uncertainties> ] alias <alias> 

Specifies the cross section of a given process <initial_state> -> 
<f inal_state_particles>, where <f inal_state_particles> is a list of particle 
names. Antiparticles are specified by <particle>~. Momentarily only e"'"e~ processes 
are implemented, identified by <initial_state>=ee. The centre-of-mass energy (if 
no unit is given, GeV is assumed) is given by <Ecms>, the polarisation of the incoming 
particles are given by <polarisationl> and <polarisation2>. If no unit of the 
cross-section and the uncertainty is given, they are assumed to be given in fb. The 
integer number <alias> identifies the cross-section in the correlationCoef f icient 
commands. No value and uncertainty must be given if the cross section is not used as 
an observable on its own but within xsbr instead. 
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Table 3: Particles known to Fittino. The following particles can be used after mass, sigma, 
BR and width. Antiparticles are identihed by ~ after the particle name. 



Particle Name 


Explanation 


Particle Name 


Explanation 


W 


W boson 


SelectronL 


selectron 


Z 


Z boson 


SelectronR 


eR selectron 


gamma 


7 


SnueL 


t'e^ electron sneutrino 


gluon 


gluon 


SmuL 


(11 smuon 


Up 


u quark 


SmuR 


jiR smuon 


Down 


d quark 


SnumuL 


v'lj,^ muon sneutrino 


Charm 


c quark 


Staul 


tl stau 


Strange 


s quark 


Stau2 


tr stau 


Top 


t quark 


Snutaul 


Pri tau sneutrino 


Bottom 


b quark 


Gluino 


g gluino 


Electron 


electron 


Neutralinol 


Xi neutralino 


Nue 


electron neutrino 


Neutralino2 


X2 neutralino 


Mu 


muon 


NeutralinoS 


X3 neutralino 


Numu 


muon neutrino 


Neutralino4 


X4 neutralino 


Tau 


tau 


Charginol 


Xi chargino 


Nutau 


tau neutrino i^r 


Chargino2 


X2 chargino 


hO 


h CP-even Higgs boson 






AO 


A CP-odd Higgs boson 






HO 


H CP-even Higgs boson 






Hplus 


H^ charged Higgs boson 






SdownL 


di squark 






SdownR 


da squark 






SupL 


ui squark 






SupR 


ur squark 






SstrangeL 


sl squark 






SstrangeR 


sr squark 






ScharmL 


cl squark 






ScharmR 


cr squark 






Sbottoml 


bi squark 






Sbottom2 


1)2 squark 






Stopl 


ti squark 






Stop2 


t2 squark 
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BR ( <decaying_particle> -> <decay_products> ) [ <value> 
+- <uncertainties> ] alias <alias> 

Specifies the branching fraction of a given particle <decaying_particle> into the decay 
products <decay_products>. The integer number <alias> identifies the branching 
fraction in the correlationCoef f icient commands. No value and uncertainty must 
be given if the branching fraction is not used as an observable on its own but within 
brsum, xsbr or brratio instead. 

brsum ( br_<aliasl> br_<alias2> [. . .] ) <value> 
+- <uncertainties> alias <alias> 

Specifies the sum of several branching fractions, which all have to be defined with a 
unique alias number before. No value and uncertainty must be given if the sum of 
branching fractions is not used as an observable on its own but within xsbr or brratio 
instead. 

width <particle> <value> +- <uncertainties> alias <alias> 

Specifies the total width of the particle <particle>. If no unit is given in <value> and 

<uncertainty>, they are assumed to be given in GeV. 

limit mass<naine> <|> <limit> 

Allows the user to specify upper or lower mass limits of yet undiscovered SUSY particles. 
The limit is used in the calculation of the of the fit in the following way: If the 
predicted mass is in agreement with the limit, the contribution of this observable to the 
total is zero. If the limit is violated by the predicted mass nip, the x^ contribution is 
((<limit>— mp)/(<limit>/10))^, i. e. a 10% violation of the limit adds 1 to the x^ of 
the fit. 

LEObs ( <naine> ) <value> +- <uncertainties> alias <alias> 

Low-energy precision observable. In Fittino 1.1 the following observables names are 

implemented: 

- bsg: BR(6 57) 

- gmin2: {g - 2)^ 

- drho: Ap 

xsbr ( signia_<aliasl> [sigma_<alias2> [...]] br_<alias3> 
[br_<alias4> [...]] ) <value> +- <uncertainties> alias <alias> 

Product of an arbitrary number of cross sections and branching fractions. This command 
must be used in the input file after all cross sections and branching fractions have been 
defined. The product of cross sections and branching fractions is then specified by 
the alias numbers of the individual observables. These are automatically not used in 
the fit anymore. Therefore the actual values entered for the individual observables are 
irrelevant, if they are used in xsbr. Instead of br_<alias> also brsum_<alias> can be 
used. 

brratio ( br_<aliasl> br_<alias2> ) <value> +- <uncertainties> alias <alias> 
Ratio of two branching fractions: br_<aliasl>/br_<alias2>. Also here the individual 
branching fractions must be defined before their use in brratio and are not used in 
the fit. Instead of br_<alias> also brsum_<alias> can be used. 

Other observables dedicated for special cases are 



14 



— sin2thetaW <value> +- <uncertainties> 

Specifies the value of sin^ 6w 

— cos2phiR <value> +- <uncertainties> 

— cos2phiL <value> +- <uncertainties> 

Specify the values of the chargino mixing angles, used for initialisation. 

Correlations among observables After all observables have been specified, the following 
command can be used to specify the correlation among observables. 

• correlationCoef f icient <observablel> <observable2> <value> 

If the observables are masses, they are identified by mass<name>. All other observables 
are identified by their alias number, such as sigma_<alias>. 

Parameters Fittino can fit any combination of the MSSM-24 and SM parameters given 
in Tab. 4 to the observables given in the input file. As alternatives to MSSM-24, high-scale 
parameters of the mSUGRA, GMSB and AMSB models can be determined. The following 
commands can be used to specify the parameters that are fitted to the observables and the 
parameters which are kept fixed. 

• fitParameter <parameter> [ <value> [ +- <uncertainty>] ] 

Specifies one of the parameters that should be fitted. The names which are avail- 
able for <parameter> are listed in Tab. 4. If a value <value> is given, then using 
FitAllDirectly it is possible to specify that <value> should be the initial value of the 
parameter in the fit. If additionally an uncertainty is given, then this can be used in the 
pull distribution calculation with CalcPullDist or CalcIndChisqContr. If no unit for 
<value> is given for a dimension-full parameter, then it is assumed to be given in GeV. 

• fixParameter <parameter> <value> 

Specifies a parameter kept fixed during the fit. 

• universality <parameterl> <parameter2> [. . .] 

Specifies that <parameter2> (and possible other following parameters) should not be 
fitted on its own, but that it should be set to the value of <parameterl> during the fit. 
This is useful if unification among generations shall be assumed. 

Flags The following flags can be used to control the behaviour of Fittino during the fit and 
to specify what operations Fittino should perform. 

• fitModel <model> 

SUSY model under study. The default value is MSSM, with the 24 low-scale SUSY 
Lagrangian parameters from Tab. 4 as parameters. If fitModel is mSUGRA, the mSUGRA 
parameters TanBeta, MO, M12 and AO can be fitted using fitParameter, and SignMu can 
be fixed with fixParameter. Thus the high-scale parameters can be fitted directly to 
the observables. Correspondingly GMSB and AMSB parameters can be fitted/fixed, if 
fitModel is set to GMSB and to AMSB respectively. In these cases the parameter names 
are TanBeta, Lambda, Mmess, cGrav, N5 and SignMu in case of GMSB and TanBeta, MO, 
M32 and SignMu for AMSB models. 
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Table 4: MSSM parameters known to Fittino. The following parameters can he used with 
fitParameter, fixParameter and universality. 



Parameter Name 


Explanation 


TanBeta 


Ratio of Higgs vacuum expectation values 


Mu 


H parameter, controls Higgsino mixing 


Xtau 


Tau mixing parameter 


Xtop 


Top mixing parameter 


Xbottom 


Bottom mixing parameter 


MSelectronR 


Right scalar electron mass parameter 


MSmuR 


Right scalar muon mass parameter 


MStauR 


Right scalar tau mass parameter 


MSelectronL 


Left 1st. gen. scalar lepton mass parameter 


MSmuL 


Left 2nd. gen. scalar lepton mass parameter 


MStauL 


Left 3rd. gen. scalar lepton mass parameter 


MSdownR 


Right scalar down mass parameter 


MSstrangeR 


Right scalar strange mass parameter 


MSbottomR 


Right scalar bottom mass parameter 


MSupR 


Right scalar up mass parameter 


MScharmR 


Right scalar charm mass parameter 


MStopR 


Right scalar top mass parameter 


MSupL 


Left 1st. gen. scalar quark mass parameter 


MScharmL 


Left 2nd. gen. scalar quark mass parameter 


MStopL 


Left 3rd. gen. scalar quark mass parameter 


Ml 


U{1)y gaugino (Bino) mass parameter 


M2 


SU{2)l gaugino (Wino) mass parameter 


M3 


SU{3)c gaugino (gluino) mass parameter 


massAO 


Pseudoscalar Higgs mass 


massW 


W boson mass 


massZ 


Z boson mass 


massTop 


Top quark mass mt(mt) 


massBottom 


Bottom quark mass mb(mb) 


massCharm 


Charm quark mass mc{mc) 
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LoopCorrections on | off 

If LoopCorrections is off, then no fit is performed but just the tree-level estimates of 
the parameters are calculated. By default it is on. 

ISR on I off 

Switches ISR corrections in the cross-section calculations on or off. By default it is on. 
UseGivenStartValues on | off 

If UseGivenStartValues is on, then the start values of the parameters in the fit are not 
determined from tree-level estimates, but from the values given in f itParameter. By 
default it is off. 

FitAllDirectly on | off 

If FitAllDirectly is on, then the initial fits of subsets of the parameter space (as 
described above in Section 2) are omitted. By default it is off. 

CalcPullDist on I off 

If CalcPullDist is on, then pull distributions for all parameters specified with 
f itParameter are calculated. It is necessary that each parameter is given with its value 
and uncertainty. The pull distribution is then calculated with respect to the parameter 
value and the width is compared with the parameter uncertainty. For each parameter, a 
ROOT [18] histogram is created in the output file PullDistributions . root. Addition- 
ally, a ROOT tree with all fitted parameters and the corresponding ^ values is created, 
which allows for the extraction of the parameter correlations. The number of fits per 
Fittino run can be specified using the command NumberPulls. By default CalcPullDist 
is off. This command is very useful to test the fitted parameters and their uncertainties 
calculated in a previous run of Fittino. 

CalcIndChisqContr on | off 

If CalcIndChisqContr is on, for each parameter specified with f itParameter the indi- 
vidual contribution of each observable to the Ax^ of the fit is calculated, if the parameter 
is varied by itlo". It is necessary that each parameter is given with its value and uncer- 
tainty. The parameter is varied once by +1 a and once by — 1 o". The resulting total A^^ 
and the individual Ax? of each observable Oj are averaged. The total Ax'^ indicates 
the correlation of the parameter with all other parameters. If the parameters are not 
too strongly correlated, the individual Ax^ provide a measure for the contribution of 
observable Oi to the determination of the parameter. The output is given in the file 
f ittino_individual_chisq_contr . out. By default, CalcIndChisqContr is off. 

BoundsOnX on | off 

If BoundsOnX is on, then the parameters X,-, Xt and X^a are bounded between Xscanlow< 
X <Xscanhigh. By default BoundsOnX is on. 

SepFitmA on | off 

If SepFitmA is on, a separate fit of m^^^^ to mA ^,1^ is performed, in order to improve 
the start value of the parameter '^Arun- This is necessary in some scenarios, where the 
difference between ttia (which is used as a tree- level estimate of rriArun) and mx^^^ is 
large and correlations of ra^^^^ are strong. By default SepFitmA is off. 

SepFitTanbMu on | off 

If SepFitTanbMu is on, a separate fit of tan (3 and /i is performed before the main fit and 
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directly after the separate fit of rriArun to mx^^^^. Only gaugino observables are used. By 
default SepFitTanbMu is off. 

ScanX on I off 

If ScanX is on, then the parameters and X\y are individually scanned in the range 
—6000 < X < 2000 before the main fit and before the separate fit of the squark sector. 
This helps to avoid local minima which typically occur at parameter values with the 
wrong sign. By default ScanX is on. 

SepFitTanbX on | off 

If SepFitTanbX is on, a separate fit of tan /?, M^^ , M^^ , X^ and Xb is performed before 
the main fit and after the separate fit of the squark sector. Only squark sector observables 
are used. By default SepFitTanbX is on. 

Calculator <calculator_name> <path> 

Specifies the tool for the calculation of the theory predictions. <path> specifies the 
location of the calculator executable. Currently SPheno is implemented, but also any 
other tool capable of input and output according to the SUSY Les Houches Accord [6] 
can be easily interfaced with Fittino. 

UseMinos on | off 

If UseMinos is on, then MINOS is used to perform a detailed uncertainty analysis after 
MINIMIZE converged. UseMinos implies automatically that UseHesse is on (see below). 
Since this can take very long (order of several days on a PHI l.SGHz in a typical fit of 
the full MSSM spectrum), this option is off by default. 

UseHesse on | off 

If UseHesse is on, then the HESSE function in MINUIT is used to perform a detailed 
error matrix calculation after MINIMIZE converged, assuming parabolic errors. By 
default, UseHesse is off. 

GetContours on | off 

If GetContours is on, after MINOS the two-dimensional uncertainty contours of 
all combinations of two parameters are calculated and stored in the output files 
FitContours . root (fit contours with an error definition of A^^ = 1 (thus the resulting 
two-dimensional surface will not be the two-dimensional 68 % region) and all parame- 
ters free) and FitContours_2params_f ree . root (just two parameters free per contour). 
Since this can take very long (order of several days to several weeks in a typical fit of 
the full MSSM spectrum), this option is off by default. 

UseSimAnnBef ore on | off 

If UseSimAnnBef ore is on, a Simulated Annealing algorithm is used directly after the 
tree-level estimates of the parameters. Typically this algorithm is very robust and finds 
the global minimum of the reliably. The maximum number of calls to the theory 
code in the Simulated Annealing process is regulated using MaxCallsSimAnn, which by 
default is 300 000. The temperature reduction can be set using TempRedSimAnn, which 
by default is 0.4. The history of the Simulated Annealing process is recorded in a ROOT 
ntuple the file SimAnnNtupFile .root. By default, UseSimAnnBef ore is off. 

UseSimAnnWhile on | off 

If UseSimAnnWhile is on and NumberOf Minimizations is larger than 1, Simulated An- 
nealing is used in between subsequent calls to MINIMIZE in the global fit. It works as 
explained for UseSimAnnBef ore. By default, UseSimAnnWhile is off. 
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NumberOf Minimizations <number> 

In very complex cases the first call to MINIMIZE often does converge near the true 
minimum of the fit, but the convergence criteria are often not fulfilled after the first call 
to MINIMIZE. Therefore, MINIMIZE can be called <number> times after each other. 
By default, NumberOf Minimizations is 1. 

ErrDef <real_number> 

If the parameter space of the fit is very complex, either due to a large number of d.o.f. 
or because of observables with large uncertainties, sometimes MINOS is not able to find 
a positive definite error matrix with the standard setting of Ax^ = 1 for the definition 
of the 1 a bound of the parameters. Therefore, using ErrDef, the error definition can be 
changed from 1 to any other positive number. After MINOS is finished, the uncertainties 
and the covariance matrix found by MINOS is re-transformed assuming parabolic errors 
such as to represent an error definition of 1. With small error definitions, MINOS finds 
the uncertainties more easily, since it is more seldom trapped in local minima close to 
the absolute minimum. On the other hand, the uncertainties on the parameters are less 
precise for small error definitions. By default, ErrDef is set to 1. 

NumberPulls <number> 

Specifies the number of individual fits for the calculation of pull distributions in one run 
of Fittino. For each fit, the observables are smeared in a Gaussian form according to 
their covariance matrix. The initialisation of the random number generator used for the 
smearing uses either the seed provided with RandomGeneratorSeed or the system time 
in seconds plus the system uptime in seconds plus the PID of the Fittino process plus 
the amount of currently available free swap space. Thus it is highly improbable that 
two Fittino processes share the same initialization. By default, NumberPulls is set to 
10. Only MINIMIZE is used in the fit, MINOS is switched off. 

TempRedSimAnn <number> 

Temperature reduction factor for the Simulated Annealing algorithm. Should be in the 
interval ]0, 1[. The default value is 0.4. 

MaxCallsSimAnn <number> 

Maximum number of calls to the theory code in one Simulated Annealing run. The 
default value is 300 000. 

InitTempSimAnn <number> 

Initial temperature of the Similated Annealing algorithm. If it is not set or negative, 
the initial temperature is calculated from the variation of the at the starting point. 

Verbose on | off 

If Verbose is on, all parameter and observable values are printed to stdout for every 
iteration of the fit. If Verbose is off, only every tenth iteration is shown. By default. 
Verbose is on. 

ScanParameters on | off 

If ScanParameters is on, either a one- or two-dimensional scan of the surface is 
performed and stored in the ROOT file ParameterScan.root. At most two parameters 
to be scanned must be specified using scanParameter. The scan is performed after the 
fit. By default, ScanParameters is off. 
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• scanParameter <paraiiieter> <lbound> <ubound> <nsteps> 

Specifies one parameter to be scanned if scanParameters is on. The lower bound 
<lbound> and upper bound <ubound> are assumed to be in GeV if no unit is given. The 
argument <nsteps> specifies the number of steps in this dimension of the scan. One or 
two scanParameter statements are needed for a scan. 

• PerformFit on | off 

If PerformFit is off, no fit is performed. Instead, evaluation tasks like scanParameters 
are executed directly. By default PerformFit is on. 

• RandomGeneratorSeed <value> 

If RandomGeneratorSeed is given, a seed <value> is used for the initialization of 
the random number generator in CalcPullDist and in the simulated annealing. If 
RandomGeneratorSeed is not given, the seed is calculated from a combination of system 
uptime, free swap space, pid of the Fittino process and system time, which should be 
unique for every run of Fittino. 

• MaxCalculatorTime <value> 

Specified the maximal time in seconds allowed for one calculation of the SUSY observ- 
ables by the calculator code. The default value is 20 seconds. 

An example for a Fittino input file fittino. in can be found in Section A.l. 

4 Quickstart Guide 

After specifying the observables with the appropriate commands from Tab. 1, the parameters 
to be fitted with f itParameter and the fixed parameters with f ixParameter, a fit with the 
following set of flags from Tab. 2 can be tried: 



LoopCorrections off 

ISR on 

UseGivenStartValues on 

FitAllDirectly on 

BoundsDnX off 

ScanX off 

SepFitTanbX off 

SepFitTanbMu off 

SepFitmA off 

UseMinos off 

UseHesse off 

GetContours off 

UseSimAnnBef ore off 

Us eSimAnn While off 

ScanParameters off 

CalcPullDist off 

CalcIndChisqContr off 
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Calculator 



YOUR_CALCULATOR /path/to/your/theory/code 



NumberOf Minimizations 1 
ErrDef 1 . 

With this parameter setting Fittino perfoms no fit but just calculates the tree-level estimates 
of the parameters, if in the Fittino code the theory program YOUR_CALCULATDR is known. 
If the code stops because a tree-level estimate could not be made (e.g. "Value massSupL 
not found"), then additionally the requested observable should be specified using the nof it 
command. This is done in order to give the user the responsibility over the inputs used for 
the tree- level estimates. The nof it command makes sure that the observable is not used in 
the fit. 

After the tree-level estimates have successfully been determined, 
LoopCorrections on 

should be set. This ensures that the fit is actually started. No Simulated Annealing or 
subsector fit is performed. If the 'true' parameter values are known, for the first try those 
should be entered in the fitParameter and f ixParameter commands, in order to ensure 
that the global fit starts with the true parameter settings. The value printed after the 
first call to the theory code should then be close to 0, if all parameters and observables 
are entered correctly. Even if the observables are smeared within their experimental and 
theoretical uncertainties, the start of the fit with the original set of parameters can be helpful 
to learn about the behaviour of the system. The reason for this step is that the start with the 
'true' parameter settings allows a fast check whether the system is constrained enough around 
the minimum to yield a realistic uncertaity estimate. If the resulting uncertainty matrix in 
fittino. out is not called Error matrix accurate, then the significance of the observables 
for the determination of the uncertainties of the parameters at the minimum is not strong 
enough. New observables have to be added, the observable uncertainties have to be reduced (if 
feasible) or the number of parameters has to be decreased. If nothing of the above is possible, 
a fit with 

UseMinos on 

may still yield a correct error matrix in many cases. 

After it is ensured that the problem is constrained strongly enough to yield a meaningful 
uncertainty matrix, the full fit can be tested. This can be done either with the subsector fit 
method or with Simulated Annealing. In the former case, 

UseGivenStartValues off 
FitAllDirectly off 

has to be set. If the subsector fits do not converge closely enough to the true minimum of 
the problem, switching on some of the additional subsector fits may help, depending on the 
individual problem. Large values of the in the early phases of the subsector fits and the 
subsequent global fit are generally no reason to worry. 

If a more stable (but also slower) way should be tried, the Simulated Annealing can be 
used, in this case 

UseGivenStartValues off 
FitAllDirectly on 
UseSimAnnBef ore on 
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has to be used. The default setting of 

TempRedS imAnn . 4 

MaxCallsSimAnn 300000 

ensures a stable convergence for most problems . If it is felt that the convergence could be 
faster, a smaller value of InitTempSimAnn than the one chosen by Fittino and maybe a smaller 
TempRedS imAnn can be tried. Should no convergence occur, all values can be increased. 

If the fit converges to a minimum close to the expected one instead of the true minimum 
(or, in case the true parameter values are unknown, it is felt that the convergence is not 
sufficient), Fittino can be made even more robust by setting NumberOf Minimizations to a 
value larger than one (typically 2 or 3), and using 

UseSimAnnWhile on 

This ensure that several subsequent global fits with intermediate Simulated Annealing algo- 
rithms are performed. 

After the correct minimum is found and the global fit successfully converged, the uncer- 
tainty matrix can be refined and the asymmetric uncertainties can be calculated using 

UseMinos on 

After the fitted parameter values and their uncertainties have been obtained in this way, a fur- 
ther analysis of the parameter determination can be done using GetContours, CalcPullDist 
or CalcIndChisqContr. In the case of CalcPullDist, the parameter uncertainties found by 
Fittino in the previous run have to be entered together with the true parameter values in 
f itParameter. In the case of CalcIndChisqContr, both the fitted central parameter values 
and the parameter uncertainties as found by Fittino have to be entered in f itParameter. 

5 The Fittino Output Files 

Depending on the requested operation, Fittino saves its results in the following files: 

• fittino. out: Main output file. It contains the observables, their covariance matrix, 
the fixed and fitted MSSM parameters and their correlation and covariance matrices. 
Additionally, information on the ^-nd the accuracy of the error matrix estimate is 
shown. An example of this output file is given in Section A. 2. 

• FitContours . root: In case a MINOS uncertainty analysis has been performed (using 
UseMinos), the two-dimensional uncertainty contours of each pair of parameters are 
calculated if the option GetContours is on. The contours are stored in the ROOT [18] 
format in this file. 

• FitContours_2params_f ree . root: As FitContours . root, but for each two dimensional 
contour of two parameters all other parameters are fixed to their best fit values. 

• f ittino_individual_chisq_contr . out: In case that an analysis of which observables 
determine which parameters is requested using CalcIndChisqContr, the individual con- 
tributions to the Ax^ are listed in this file. 

• PullDistributions . root: In case pull distributions are calculated using 
CalcPullDistr, the pull distributions and the distribution of the total of tlis 
fits are stored in this file in the ROOT format. One-dimensional projections of the 
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distribution of fitted parameter values are available in ROOT histograms. The full 
correlation of parameters can be deduced from the ROOT tree containing all fitted 
parameter values and the corresponding values. 

• SimAnnNtupFile . root: This file contains one ROOT-ntuple for each run of the Simu- 
lated Annealing algorithm. The ntuple contains the value and the parameter settings 
at each step. 

• ParameterScan.root: In case a parameter scan is performed, a TGraph (or TGraph2D) 
object describing the curve (plane) is stored in this file. 

Examples for the contents of the output files are given in the appendix. 

6 Summary 

Fittino is a versatile program with flexible user interface which is designed to perform a 
global fit of MSSM and SM parameters to observables from present and future collider 
experiments. Using Fittino, information about the feasibility of a unique reconstruction of 
MSSM parameters can be obtained, parameter uncertainties and correlations for simulated sets 
of observables can be determined and it can be used to optimise the observables for a desired 
parameter uncertainty. Additionally, the fit results of Fittino can be tested by the automatic 
generation of pull distributions and distributions. Two-dimensional uncertainty contours 
can be obtained to study nonlinearities in the parameter correlations and to visualise the 
results. The contribution of the individual observables to the confinement of each parameter 
can be obtained by studying the individual contributions of each observable to the Ax^ for 
parameter variations by iblcr. 

Fittino has been tested using global fits of 19 MSSM parameters and one SM parameter to 
observables from the LHC and the ILC. The program has also been used to directly fit high- 
scale mSUGRA parameters to anticipated LHC and ILC measurements. No prior knowledge of 
the parameters has been assumed at any step. The parameters are reconstructed entirely from 
the observables. It has been shown that the fit strategy of Fittino, using tree-level estimates 
of the parameters, subsector fits in parts of the parameter space to improve the tree- level 
estimates and then global fits is stable even for different MSSM scenarios like SPSla, SPSla' 
and SPS7. The possibility to include Simulated Annealing in order to find the global x^ 
minimum before the global fit provides a further stable method to find the optimal parameter 
values. 

In the future, Fittino can be used for systematic studies of the dependence on measurement 
uncertainties. The most crucial observables for the parameter determination can be identified 
and the detector and machine, especially the distribution of the luminosities with different 
centre-of-mass energies and beam polarisations, can be optimised for most precise parameter 
measurements. Also different sets of SUSY parameters and their phenomenology can be 
studied. The infiuence of theoretical uncertainties can be determined and regions with need 
for improvement can be identified. 
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A Fittino Input and Output Files 
A.l Example for the Fittino Steering File 

In the following an example for the Fittino input file fittino . in is given. For demonstration, 
it shows a fit of just two observables, tan (3 and fj,. All other parameters are kept fixed. An 
assorted selection of input observables from the ILC is used for illustration. The order of the 
inputs is mostly arbitrary, apart from the correlations of the observables, which always must 
be stated after all observables are given. 

######################################################################### 

### Fittino input file ### 

######################################################################### 

### This is an example steering file. ### 

### P. Bechtle, 20040520 ### 

######################################################################### 



# SM observables 

massW 80.3078 GeV +- 0.039 GeV 

massZ 91.1187 GeV +- 0.0021 GeV 

massTop 174.3 GeV +- 0.05 GeV 

massBottom 4.2 GeV +- 0.5 GeV 

sin2ThetaW 0.23113 +- 0.00015 

alphas 0.1172 +- 0.002 



# Higgs sector at ILC 500 and ILC 1000 

masshO 108.67 GeV +- 0.05 GeV +- 0.5 GeV # Degrassi et al 



massAO 


311 


.036 


GeV 


+- 


1.3 GeV 




# not : 


# Gauginos at ILC 500 


and ILC 1000 










massNeutralinol 


95.9804 


GeV 


+- 


0.05 GeV 


# 


~chi_10 


massNeutralino2 


180 


.397 


GeV 


+- 


0.08 GeV 


# 


~chi_20 


massCharginol 


180 


.033 


GeV 


+- 


0.55 GeV 


# 


~chi_l+ 


massChargino2 


381 


.929 


GeV 


+- 


3.0 GeV 


# 


~chi_2+ 


# plus squark and slepton masses 


, if 


necessary with preceed: 


nofit massSupL 


561 


.539 


GeV 


+- 


9.8 GeV 


# 


~uL 


nofit massSupR 


543 


.35 


GeV 


+- 


23.6 GeV 


# 


~uR 


nofit massSdownL 


561 


.54 


GeV 


+- 


9.8 GeV 


# 


~dL 


nofit massSdownR 


543 


.348 


GeV 


+- 


23.6 GeV 


# 


~dR 


nofit massSstrangeL 


561 


.54 


GeV 


+- 


9.8 GeV 


# 


~sL 


nofit massSstrangeR 


543 


.348 


GeV 


+- 


23.6 GeV 


# 


~sR 


nofit massScharmL 


561 


.54 


GeV 


+- 


9.8 GeV 


# 


~cL 


nofit massScharmR 


543 


.348 


GeV 


+- 


23.6 GeV 


# 


~cR 


nofit massSbottoml 


502 


.059 


GeV 


+- 


5.7 GeV 


# 


~b_l 


nofit massSbottom2 


541 


.81 


GeV 


+- 


6.2 GeV 


# 


~b_2 


nofit massStopl 


365 


.819 


GeV 


+- 


2.0 GeV 


# 


~t_l 


nofit massStop2 


600 


.037 


GeV 


+- 


20.0 GeV 


# 


~t_2 


nofit massSelectronL 


190 


.209 


GeV 


+- 


0.2 GeV 


# 


~e_L- 


nofit massSelectronR 


124 


.883 


GeV 


+- 


0.05 GeV 


# 


~e_R- 


nofit massSnueL 


172 


.947 


GeV 


+- 


0.7 GeV 


# 


~nu_eL 


nofit massSmuL 


190 


.237 


GeV 


+- 


0.5 GeV 


# 


~mu_L- 
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nof it 


massSmuR 


124. 


837 


GeV 


+- 


0. 


2 GeV 


# 


~inu_l 


R- 


nof it 


massStaul 


107. 


292 


GeV 


+- 


0. 


3 GeV 


# 


~tau. 


_1- 


nof it 


massStau.2 


195. 


290 


GeV 


+- 


1. 


1 GeV 


# 


~tau. 


_2- 


nof it 


massGluino 


603. 


639 


GeV 


+- 


6. 


4 GeV 


# 


o 




nof it 


massNeutralinol 


97.7662 


GeV 


+- 


0. 


05 GeV 


# 


~chi. 


_10 


nof it 


massNeutralino2 


184. 


345 


GeV 


+- 


0. 


08 GeV 


# 


~chi. 


_20 


nof it 


massNeutralino3 


-404.134 GeV +■ 


- 4 


:.0 GeV 


# 


~chi. 


_30 


nof it 


massNeutralino4 


417. 


037 


GeV 


+- 


2. 


3 GeV 


# 


~chi. 


_40 


nof it 


massCharginol 


184. 


132 


GeV 


+- 


0. 


55 GeV 


# 


~chi. 


_1+ 


nof it 


massChargino2 


418. 


495 


GeV 


+- 


3. 


GeV 


# 


~chi. 


_2+ 



# possible edges 

# edge 1 massNeutralinol massNeutralino2 263.50279 GeV +- 1.2 GeV alias 1 

# edge 2 massNeutralinol massNeutralino2 79.09719 GeV +- 1.2 GeV alias 2 

# estimated chargino mixing angles from ILC 500 

nof it cos2PhiL 0.6737 +- 0.05 # rough estimate 

nof it cos2PhiR 0.8978 +- 0.05 # rough estimate 

# ILC 500 Cross-sections 

Sigma ( ee->Neutralinol Neutralino2 , 500 . , . 8 , -0 . 6 ) 20.5026 fb +- 2.0 fb alias 1 

Sigma ( ee->Neutralino2 Neutralino2 , 500 . , . 8 , -0 . 6 ) 5.62767 fb +- 2.0 fb alias 2 

Sigma ( ee->Charginol Charginol~, 500. ,0.8,-0.6 ) 13.6598 fb +- 1.0 fb alias 6 

Sigma ( ee->Z h0,500., 0.8, -0.6 ) 29.0566 fb +- 0.21 fb alias 7 

Sigma ( ee->Charginol Charginol~, 500. ,-0.8,0.6 ) 462.321 fb +- 5.0 fb alias 8 

Sigma ( ee->Neutralinol Neutralino2 , 500 . , -0 . 8 , . 6 ) 180.84 fb +- 2.0 fb alias 9 

Sigma ( ee->Neutralino2 Neutralino2 , 500 . , -0 . 8 , . 6 ) 200.396 fb +- 2.0 fb alias 10 



# ILC 500 Branching Fractions 

# BR ( hO -> Bottom Bottom" ) 

# BR ( hO -> Charm Charm" ) 

# BR ( hO -> Tau Tau" ) 



0.824057 +- 0.019 alias 1 

0.0405547 +- 0.02 alias 2 

0.134444 +- 0.02 alias 3 



# Product of cross sections and branching fractions: 

# xsbr ( sigma_7 br_l ) 32.874 +- 0.5 alias 1 

# Ratio of branching fractions 

# brratio ( br_l br_2 ) 20.6259 +- 0.5 alias 1 



# Correlations among observables 

# correlationCoef f icient massCharginol massNeutralinol 0.05 



# Parameters to be fitted 

fitParameter TanBeta 10.0 

f itParameter Mu 358.635 GeV 



# Fixed Parameters 

fixParameter Atau -3884.46 GeV 

fixParameter MSelectronR 135.762 GeV 

fixParameter MStauR 133.564 GeV 



25 



f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 
f ixParameter 



MSelectronL 

MStauL 

Atop 

Abottom 

MSdownR 

MSbottomR 

MSupR 

MStopR 

MSupL 

MStopL 

Ml 

M2 

M3 

mass AO 
massTop 
massBottom 
mass Charm 



195.21 GeV 
194.277 GeV 
-506.936 GeV 
-4444.56 GeV 
528.135 GeV 
524.719 GeV 
530.244 GeV 
424.515 GeV 
548.704 GeV 
499.986 GeV 
101.814 GeV 
191.771 GeV 
588.798 GeV 
399.763 GeV 
174.3 GeV 
4.200e+00 GeV 
1.2e+00 GeV 



# Universalities among parameters 
f ixParameter MSmuR 
f ixParameter MSmuL 
f ixParameter MSstrangeR 
f ixParameter MScharmR 
f ixParameter MScharmL 



135.76 GeV # 1.35760e+02 GeV # 135.76 GeV 
195.21 GeV # 1.95199e+02 GeV # 195.21 GeV 
528.14 GeV # 5.28134e+02 GeV # 528.14 GeV 
530.253 GeV # 5.30270e+02 GeV # 530.253 GeV 
548.705 GeV # 5.48702e+02 GeV # 548.705 GeV 



# possible unifications 

# universality MSelectronR MSmuR 

# universality MSelectronL MSmuL 

# universality MSdownR MSstrangeR 

# universality MSupR MScharmR 

# universality MSupL MScharmL 



# switches 

LoopCorrections on 
ISR on 



# Use full loop corrections in SPHENO 

# Switch on ISR 



UseGivenStartValues on 



FitAllDirectly 



on 



ScanParameters off 
CalcPullDist off 
CalcIndChisqContr off 



# Calculator 



SPHENO /home/bechtle/programs/SPheno2 . 2 . 2/SPheno 



UseMinos 



on 



NumberOf Minimizations 1 
ErrDef 1 , 
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NumberPulls 



3 



A. 2 Example for the Fittino Output File 

The result of the fit from the input file fittino. in shown in Section A.l is given in the file 
fittino. out. It is presented in the following. The lengthy table of the correlations of the 
observables is abbreviated. Please note the message Error Matrix Accurate in the end of 
the file. If this message is not given, the fit uncertainties are not reliable, since MINUIT could 
not find a positive definite uncertainty matrix. 

#################### Fittino Fit Summary #################### 

created by Fittino version 1.0.3 

on Wednesday, August 10, 2005 at 07:28:02 

Input values : 



ee 
ee 
ee 

ee 
ee 
ee 



massZ 
massW 
massCharm 
massBottom 
massTop 
massTau 
alphas 
alphaem 
sin2ThetaW 
masshO 
massAO 
massNeutralinol 
massNeutralino2 
massCharginol 
massChargino2 
-> Neutralinol Neutral 
-> Neutralino2 Neutral 
-> Charginol Charginol 

ee -> Z hO 
-> Charginol Charginol 
-> Neutralinol Neutral 
-> Neutralino2 Neutral 



91.1187 
80.3078 
1.2 
4.2 
174.3 
1.77699 
0.1172 
127.934 
0.23113 + 
108.67 + 
311.036 + 
95.9804 + 
180.397 + 
180.033 + 
381.929 + 
ino2 
ino2 



0.0021 
0.039 
0.2 
0.5 
0.05 
0.00029 
0.002 
0.027 
0.00015 
0.502494 
1.3 
0.05 
0.08 
0.55 
3 

20.502600 +- 2.000000 
5.627670 +- 2.000000 



13.659800 +- 1.000000 
29.0566 +- 0.21 

462.321000 +- 5.000000 
ino2 180.840000 +- 2.000000 

ino2 200.396000 +- 2.000000 



Covariance matrix for input value: 



masshO 
mass AO 

all observables 



masshO 
0.252500 
0.000000 



massAO 
0.000000 
1.690000 



all observables 
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Fixed values: 



Xtau 
MSelectronR 
MStauR 
MSelectronL 
MStauL 
Xtop 
Xbottom 
MSdownR 
MSbottomR 
MSupR 
MStopR 
MSupL 
MStopL 
Ml 
M2 
M3 

massAO 
MSmuR 
MSmuL 
MSstrangeR 
MScharmR 
MScharmL 
massTop 
massBottom 
massCharm 



-3884.46 
135.762 
133.564 
195.21 
194.277 
-506.936 
-4444.56 
528.135 
524.719 
530.244 
424.515 
548 . 704 
499.986 
101.814 
191.771 
588.798 
399.763 
135.76 
195.21 
528.14 
530.253 
548.705 
174.3 
4.2 
1.2 



Fitted values: 



TanBeta 
Mu 



9.99992 +- 
358.639 +- 



0.0935112 
1 . 00434 



Covariance matrix for fitted parameters: 



TanBeta Mu 
TanBeta 0.00874435 -0.0532778 
Mu -0.0532778 1.00871 

Correlation matrix for fitted parameters: 



TanBeta Mu 
TanBeta 1 -0.567283 

Mu -0.567283 1 



Chisq of the fit ; 
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Figure 4: Examples of the two-dimensional uncertainty contours of the Fittino SPSla' fit. 



chisq = 0.000007 
Status of the minimization: 



Error Matrix accurate 
##################### End of Fit Summary #################### 

A. 3 Examples for the Two-Dimensional Uncertainty Contour Plots 

If the flag UseMinos is set, two sets of two-dimensional uncertainty contours of all possible 
combinations of all parameters are produced. In the first set, contained as two-dimensional 
TGraph objects in the output file FitContours_2params_f ree . root, all parameters but the 
two plotted ones are kept fixed. In the second set, contained in FitContours .root, all pa- 
rameters are free, yielding larger allowed areas inside the plotted Ax^ = 1 area. Fig. 4 shows 
examples of the obtained curves for a fit to the benchmark scenario SPSla, which are ellipses 
in the case of almost linear dependences of the observables on the parameters near the fitted 
minimum. The slope of the ellipses indicates the parameter correlations. 



A. 4 Examples for the Pull Distributions 

If parameter values and uncertainties (as determined with Fittino) are supplied for each pa- 
rameter specified with f itParameter, pull distributions and the distribution of the fit is 
automatically calculated by Fittino if the flag CalcPullDist is set. Examples for these dis- 
tributions, as obtained in a fit in the SPSla scenario with 20 parameters, can be found in 
Figures 5 and 6. 
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Figure 5: Pull distributions of the parameter tan/3 of the SPSla' fit for 120 independent fits 
with observables smeared within their uncertainties. The uncertainties of all parameters are 
well described. 




Figure 6: Distribution of for the SPSla' fit. 273 fits with 71 d .o.f. are performed. The 
observables of the simulated measurements are smeared in each fit according to their uncer- 
tainties. The mean of the distribution is Xmean — 69.2, in agreement with the expectation of 
71 ±0.9. 
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A. 5 Examples for the Observable Importance Determination 

In the observable importance determination (using CalcIndChisqContr), the individual con- 
tribution of each observable to the determination of each parameter is calculated. Addi- 
tionally, the total Ax^ with respect to the minimal for a variation of it lex of each pa- 
rameter is shown. If the total A^^ is large, then the parameter is strongly correlated with 
other parameters. If it is close to 1, then there is hardly any correlation. In the output file 
f ittino_individual_chisq_contr . out, for each parameter first all individual Ax^ values 
from each observable are shown, followed by the percentage of the contribution from the 5 
largest contributors. For the cross sections, the alias number is also shown in order to ensure 
the identification of cross sections at different polarisations and centre-of-mass energies. 



Individual Delta Chisq contributions to Parameter TanBeta (10.0005 +- 0.333429) 



list of all observables 



Delta chi~2 



massNeutralinol 
massNeutralino2 
massCharginol 
massChargino2 
ee -> Neutralinol Neutralino2 (19) 
ee -> Neutralino2 Neutralino2 (20) 
ee -> SelectronL SelectronL~ (23) 
ee -> SmuL SrnuL" (25) 
ee -> Staul Staul" (27) 
ee -> Staul Staul~ (28) 
ee -> Charginol Charginol (41) 
ee -> Z hO (45) 
ee -> Z hO (46) 



0.843007 
2.70032 
0.0617443 
0.000608651 
0.0101081 
0.0636789 
0.00277659 
0.000133514 
0.000536534 
0.00530975 
0.0489321 
0.00190326 
0.0012801 



percentage of strongest contributions: 



massNeutralino2 
ee -> Neutralino2 Neutralino2 (20) 
ee -> Charginol Charginol (41) 
massNeutralinol 
ee -> Neutralinol Neutralino2 (19) 
total Delta Chisq 



0.280837 
0.191446 
0.126957 
0.0876738 
0.0717215 
9.615270 



Individual Delta Chisq contributions to Parameter Mu (358.644 +- 1.14432) 
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